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ABSTRACT 

We made stereoscopic observations of the Vela Pulsar region with two of the 10 m diameter 
CANGAROO-III imaging atmospheric Cherenkov telescopes in January and February, 2004, in a 
search for sub-TeV gamma-rays from the pulsar and surrounding regions. We describe the observa- 
tions, provide a detailed account of the calibration methods, and introduce the improved and bias-free 
analysis techniques employed for CANGAROO-III data. No evidence of gamma-ray emission is found 
from either the pulsar position or the previously reported position offset by 0.13 degree, and the 
resulting upper limits are a factor of five less than the previously reported flux from observations 
with the CANGAROO-I 3.8 m telescope. Following the recent report by the H.E.S.S. group of TeV 
gamma-ray emission from the Pulsar Wind Nebula, which is ~0.5 degree south of the pulsar position, 
we examined this region and found supporting evidence for emission extended over ~0.6 degree. 
Subject headings: gamma rays: observations — pulsars: individual (Vela) — methods: data analysis 



1. INTRODUCTION 

The Vela Pulsar is t he brightest object in the 
sky at 100 MeV energies ijHartman et ; all Il999fl, with 
emiss ion extending beyond" 10 GeV ijTriomDSon et all 
120051). The emission at these energies is totally pulsed 
l)Kanbach et alJ 119941). The pu lsar undergoes large 
glitches (e.s.. iDodson et aTll2002|) . and has one of the 
highest values of E/d 2 of all pulsars, making it a prime 
target for s outhern hemisph ere searches for TeV gamma- 
rays (e.g., iNel et all Il992[) . Early searches for TeV 
gamma-ray emission relied on the detection of a pulsed 
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signal, and produc ed a number of sugg estive but incon- 
sistent results (see Edward s et al.lll994L for a review) . 

The first search sensitive to both pulsed and steady 
emission and made using the proven imaging atmo- 
spheric Cherenkov technique was undertaken with the 
CANGAROO-I 3.8 m telescope. A total of 119 hours 
of usable on-source data taken between 1993 and 1995 
yielded a 5.8cr excess, corresponding to a flux of 
(2.9±0.5±0.4) x 10~ 12 cm^s- 1 above 2.5±1.0TeV, 
offset from the Vela Pulsar by 0.13 degrees and 
with no significant m odulation at the pulsar period 
(Yoshikos hi et all 11997). A further 29 hours data from 
CANGAROO-I observations in 1997 showed a 4.1ct ex- 
cess at a posi tion consistent with peak from the earlie r 
observations (j Yoshikoshi et al.lll997t IYoshikoshilll998(K 
Observations by the Durham group in 1996 produced 
3cr upper limits above 300 GeV of 5 x 10~ n and 1.3 x 
10 _11 cm _2 s _1 for steady an d pulsed emission, re spec- 
tively, at the pulsar position ijChadwick et al.l|2000|) . No 
significant excess was found in a 2°x2° area centered at 
the pulsar position. 

Based on independent M onte-Carlo simulations of th e 
CANGAROO - I tele scope i|Dazelev fc Patterson! l20"o"l|l 
iDazelev et all l)2001|) searched for a TeV gamma-ray sig- 
nal in the 29 hour 1997 dataset , but found no evidence of 
emission. Da zelev et all l|2001fl noted that increasing the 
value of the length cut toward the value used in the ear- 
lier analysis did result in a somewhat increased excess, 
though not to a significant level, and acknowledged that 
the simulations had not included several known charac- 
teristics of the CANGAROO-I telescope hardware. An 
analysis of data from the Crab Nebula taken between De- 
cember 1996 and January 1998 with two different sets of 
cuts a lso failed to yield a significant excess l)Dazelev et alJ 
1200 1|) . although the resulting upper limit was higher 
than the established flux of the Crab Nebula. The TeV 



2 



Enomoto et al. 



flux from the Crab N ebula from earlier CANGAROO-I 
ijTanimori et alll998H obser vations is consistent with the 
most recent measurements (Eh aroman eTaIll2004bl) . 

The first preliminary results from observations with 
the H.E.S.S. stereoscopic imaging Cherenkov telescopes 
showed no evidence of emission from the Vela Pulsar 
nor fr om previous CANGAROO-I position llKhelifi et alJ 
|2005|) . However, very recently the H.E.S.S. group re- 
ported the detection of broad TeV emission ^0.5 degree 
south from the pulsar position e xtended over a region < 
0.6 degree (Khclifi et al. 2005b). The preliminary non- 
detection of emission from the pulsar was confirmed, yet 
the extended emission from the region of the Vela pul- 
sar win d nebula was at the level of 0.5 crab (Khclif Tet al.l 
2005b) . Here we report the results from our observations 
of the Vela Pulsar region made in early of 2004 with the 
new CANGAROO-III stereoscopic imaging atmospheric 
system. 

The previous CANGAROO-I analyses used slightly dif- 
ferent cut values, especially in the image analyses, due in 
part to differences in source declination and sky bright- 
ness. Analysis using three or four image parameters re- 
quires selection of six to eight cut parameters (i.e., both 
upper- and lower-cut values for each parameter). While 
the cuts were guided by simulations, which endeavored 
to incorporate the charact erstics of all telescope hard - 
ware and electronics (e.g., iDazelev fc Patterson! 120011. 
any fine tuning introduces a number of degrees of free- 
dom into the analysis, which can be difficult to account 
for a posteriori. The number of parameters can be re- 
duced to one w ith the use of mathem atical methods such 
as Likeliho od llEnomoto et aP 120021 or the Fisher Dis- 
criminant (Fisher 1936). In this report, we emphasize 
the application of such methods in order to remove any 
potentially bia sed operation from the analysis (see. e.g., 
iWeekes ll2005j) . Previous CANGAROO-II analyses have 
been carried out using the Likelih ood ratio, a Likelihood 
analysis with only one image cut (jEnomoto et al.l l20021 . 
The Fisher Discriminant analysis introduced here now 
removes all cut uncertainties from the image analysis. 

In §2 we describe the CANGAROO-III system, and 
in §3 we present the various calibration checks made to 
confirm its performance, intr oducing unbiased ana lysis 
methods based on Likeli hood llEnomoto et alJl2002H and 
the Fisher Discriminant l)Fisherlll936|1 . In SS4 and 5 the 
results of the CANGAROO-III observations of the Vela 
Pulsar and the surrounding field are presented, with con- 
clusions given in §6. 

2. THE CANGAROO-III STEREOSCOPIC SYSTEM 

The imaging atmospheric Cherenkov technique 
(IACT) was pioneered with the Whipple group's 
detec tion of TeV gamm a-rays from the Crab Neb- 
ula ijWeekes et alJ Il989|) . This technique enabled 
TeV gamma-rays to be distinguished from the huge 
background of cosmic rays with t he use of th e "Hillas 
moments" of the Cherenkov images llHillasll985l1 . Stereo- 
scopic observations were successfully demonstrated by 
the HEGRA group to significan tly improve the rejectio n 
of the cosmic ray background (jAharonian et al.l Tl9991. 
More recently, the H.E.S.S. group has reported the 
detection of a number of faint gamma-ray sources 
with angular resolution s of only a few-arc minutes 
ijAharonian et alJl2005a^ . 



The CANGAROO-III stereoscopic system consists of 
four imaging atmospheric imaging telescopes located 
near Woomera, South Australia (31°06'S, 136°47'E, 160 
m a.s.L). Each telescope has a 10m diameter reflec- 
tor consisting of 114 spherical mirror segments, each of 
80 cm diameter and with an average radius of curvature 
of 16.4 m. Th e segments are made of Fiber Reinforced 
Plastic (FRP) (jKawachi et all2001|1 and are aligned on a 
parabolic frame with a focal length of 8 m. The total light 
collecting area is 57 m 2 . The four telescopes are located 
on the corners of a diamond with sp acings of approx- 
imately 100 m l|Enomoto et alJl2002bH . The first tele- 
scope, designated Tl, is the CANGAROO-II telescope, 
and data from this telescope was not used in the analysis 
presented here as its smaller field of view (FOV) was less 
well suited for these stereoscopic observations. The sec- 
ond and third CANGAROO-III telescopes, T2 and T3, 
were used for the observations described here, with the 
fourth telescope (T4) becoming operational after these 
observations were completed. The camera systems for 
T 2, T3, and T4 are ide ntical and are described in detail 
in iKabuki et alJ ((2003) . The pixel timing resolution is 
Ins at 20 photoelectron (p.e.) input. The noise level is 
0.1 p.e. The relative gains of pixels are adjusted to within 
10% at the hardware level by adjusting the high voltages 
applied to photomultipliers. The offline calibration are 
carried out using LED flasher data which were taken once 
per day. The gain uniformity of pixels is believed to be 
less than 5%. 

The observations were made using the so-called "wob- 
ble mode" in which the pointing position of each tele- 
scope was alternated between ±0.5 degree in declina- 
tion from the cen ter of the target every 20 minutes 
l|Daum et alJ 1199711 . Two telescopes, T2 and T3, were 
operational at this time, with T3 located 100 m to the 
south-south-east of T2. The data from each telescope 
were recorded with the trigger condition requiring that 
more than four pixels exceeded 7.6 p.e. The GPS time 
stamp was also recorded for each event. In the offline 
analysis, we combine these datasets by comparing the 
GPS timing. The distribution of timing differences is 
shown in Fig. ^ The broadness of the distribution was 
due to the time resolution of the GPS clocks. We re- 
quired an off-line coincidence between the two telescopes' 
trigger times within ±100 /is. Events outside this time 
window were considered to be accidental coincidences. 
The trigger rates of the individual telescopes were at 
most 80 Hz and this was reduced to 8 Hz by requiring 
the above coincidence. 

3. CALIBRATION 

In order to determine the sensitivity of the telescope, 
the efficiency of the light collecting system is the most 
important parameter to be precisely measured. Electro- 
magnetic showers can be accurately simulated by com- 
puter code and cosmic ray showers can be observed in 
OFF-source observations and generated using Monte- 
Carlo codes such as GEANT. We measured the point 
spread function (PSF) and reflectivity of each mirror seg- 
ment during their production and the quantum efficiency 
of each photomultiplier (PMT) pixel was also measured. 
Some deterioration in these values after the construction 
of the telescope is inevitable, and was monitored by reg- 
ular PSF measurements of bright stars. Light collection 
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Fig. 1. — GPS time differences in units of /is for all events. 



efficiencies were monitored via optical systems such as 
photo-diodes and commercial CCD cameras observing 
stars both directly, and after reflection by the mirrors. 
These showed gradual degradation in the form of broad- 
ening spot sizes and decreasing reflectivities. However, 
these data only relate to the reflectors, and we need to 
know the total performance of the whole system includ- 
ing the cameras. 

In order to measure these, we need "standard candles" ; 
sources of Cherenkov photons or non- variable gamma-ray 
sources such as the Crab Nebula. Here, we introduce our 
calibration procedure using these methods. 

3.1. Muon-ring Analysis 

The ideal light source is one that produces the same 
Cherenkov radiation that the telescope is designed to de- 
tect. Sub- TeV cosmic rays produce hadronic showers in 
the upper atmosphere at a depth of around three inelas- 
tic interaction lengths, corresponding to an altitude of 
^10 km. Secondary hadrons are produced and interact 
further, with only gamma-rays, electrons and muons sur- 
viving to near ground level. Of these, relativistic muons 
radiate Cherenkov rings near the Earth's surface, and 
these produce distinct complete or partial ring images 
on the camera. 

The contribution of these rings to the trigger rate of our 
telescopes is critical, because small Cherenkov images re- 
semble gamma-ray images. Local muon ring images can, 
however, be removed by requiring the coincidence of two 
or more telescopes separated by ~100m. Here, for cali- 
bration, we select images produced by muons at altitudes 
of 100-200 m. This can be done by selecting a somewhat 
larger "arc-length" . Arc-length is measured after fitting 
a circle to the images, and is the fitted radius multiplied 
by the recorded opening angle of the image. If larger 
images are required, we can simply restrict the selection 
to geometrically closer images such as those with alti- 
tudes < 100-200 m. The inverse of the radius (curvature) 
should have a Gaussian distribution, so we therefore plot 
the curvature as a solid histogram in Fig. [2J A peak 



Fig. 2. — Curvature distribution for the circular fitting results 
of all events for T3. The solid histogram is for all events and the 
dashed one for the selected events (with the vertical scale multiplied 
by 5). 
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Fig. 3. — Curvature distribution for the circular fitting results 
on all events for T2. The solid histogram is for all events and the 
dashed one for the selected events (with the vertical scale multiplied 
by 5). 



can be seen around 1/1.3 [1/degree] which corresponds 
to the inverse of the Cherenkov angle initiated by rela- 
tivistic muons at normal temperature and pressure. The 
distribution for T2 (which was constructed before T3), 
however, was not as good, indicating some degradation 
in its performance, as shown by the solid histogram in 
Fig. |21 These data were taken in 2003 December. We 
need to know our light collection efficiency in each ob- 
serving period, i.e., we need to improve the S/N ratio to 
select a pure muon-ring sample only. The following are 
the cut criteria used to achieve this. 

• Hit threshold more than 1 photo-electron (p.e.). 

• TDC range within ±30 ns from the mean arrival 
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time of the event. 

• At least 15 triggered (or "hit") PMTs, with a 
clustering requirement that there be at most two 
non-triggered PMTs between nearest-neighbor hit 
PMTs. 

• A circle, or ring segment, can be fitted to the image. 

• The circle has an "arc-length" of more than 2 de- 
grees. 

• A x 2 per hit-pixel, normalized to the pixel size, of 
less than 1.5. 

The curvature distributions of the samples meeting these 
criteria are shown as dashed histograms in Figs. and 
131 Muon ring events are selected with a good S/N ra- 
tio. Geometrically, the arc-length cut restricts selection 
to events occurring at < 100-200 m. The statistics of the 
accepted events are sufficient for calibration to be under- 
taken with only several hours data. 

There is some dependence of light yield on the at- 
mospheric temperature; however, this is monitored and 
recorded every second. The night-time temperature near 
Woomera ranges from near C in winter to over 30 C in 
summer. The systematic error, which includes the reflec- 
tivity of the present mirror, the PSF, and the quantum 
efficiency of the PMTs, is thought to be at the 5% level 
of the total light collection efficiency. The light collection 
efficiencies can be derived from the "size/arc-length" dis- 
tribution and the PSFs by the x 2 °f the ring fits. The 
observation periods for the Crab Nebula and Vela Pulsar 
data considered here were between the end of 2003 and 
the beginning of 2004. 

The calibrated light collection efficiencies for T2 and 
T3 were both 70%, with 5% systematic errors. These 
values are the ratios of the efficiencies derived during 
the observation period and those measured at the pro- 
duction time, i.e., the deterioration factors. The PSFs 
were approximated by Gaussian distributions with stan- 
dard deviations of 0.14, and 0.12 degrees, respectively, 
which are somewhat larger than the results obtained via 
bright star measurements at the initial installation time. 
(The corresponding values for T4, which was just com- 
ing on-line, were 85% and 0.09 degrees, respectively.) 
The x 2 distributions for the experiments (the solid his- 
togram) and the tuned Monte-Carlo simulations (the 
hatched area) are shown in Fig. 0] for T2 (upper panel) 
and T3 (lower panel). The PSF of T2 is worse than that 
of T3. The size/arc- length distributions for the experi- 
ments (the solid histogram) and the tuned Monte-Carlo 
(the hatched area) are shown in Fig. [3] for T2 (upper 
panel) and T3 (lower panel). The PSFs obtained for 
the three telescopes are somewhat larger than those of 
H.E.S.S., corresponding to a relatively lower cosmic-ray 
rejection efficiency. The resulting S/N ratio for gamma- 
rays is discussed in the next section; however, it is suf- 
ficient to detect gamma-ray sources with flux levels flux 
of ^0.1 Crab in several tens of hours. 

The measured efficiencies do not depend strongly on 
climate or atmospheric conditions, as they are local 
(near-surface) phenomena. Thus uncertainty due to Mie 
scattering, for instance, still remains. This is thought to 
be significant at the 10% level; however, the average ef- 
fect can be gauged in the Crab Nebula data described in 
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Fig. 4. — Distributions of \ 2 f° r the experimental data and the 
tuned Monte-Carlo. The solid histograms are the experimental 
data and the hatched the Monte-Carlo. The upper panel is for T2 
and the lower T3. 
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Fig. 5. — Distributions of size/arc-length for the experimental 
data and the tuned Monte-Carlo. The solid histograms are the 
experiment and the hatched the Monte-Carlo. The upper panel is 
for T2 and the lower T3. 



the following section. In fact, data taken during cloudy 
nights showed a better S/N ratio in the curvature distri- 
bution. As the muon rings are a local phenomenon, there 
is naturally no dependence on the elevation angle of the 
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Fig. 7. — Definitions of parameter Wt and Dt which are used 
in the fitting procedure described in the text. 



Fig. 6. — 9 2 distributions from Monte-Carlo simulations. The 
dashed line is obtained without the fitting procedure described in 
the text and the solid line is that with it. 



observation, whereas the cosmic ray rate varies strongly 
with elevation. Using these characteristics, we were able 
to measure the degradation as a function of time for the 
three telescopes (T2, T3, and T4), with the light collec- 
tion efficiency being found to decrease by 5% per year. 
Hereafter in the analyses of the individual sources, we 
used muon-ring data from the corresponding periods to 
tune our Monte-Carlo code. 

3.2. Crab Analysis 

The description on the Crab data used in the following 
section can be found in Appendix 1X1 At first, the Hillas 
moments were calculated for each telescope's image. The 
intersection of the major axes of the two images is the 
incident direction of the gamma-ray in the stereoscopic 
observation. The 9 2 distribution of the Monte-Carlo data 
for this observation condition calculated from the inter- 
section points is plotted in Fig. H3 (the dashed line). The 
angular resolution is, as expected, several times worse 
than that of H.E.S.S. ijAharonian et alJl2Q04T) . due to the 
larger PSFs described above. The Crab Nebula observa- 
tions were carried out at relatively low elevation angles 
and the opening angle between the T2 and T3 images are 
typically small, resulting in many parallel images. This 
results in an increased uncertainty in the intersection 
point of the two images. To avoid this, we reanalyzed the 
data with a constraint on the distances between images' 
center of gravity and the intersection point. (H.E.S.S. 
observations are made at elevation angles typically 10 
degrees higher and with smaller mirror PSF, resulting a 
good S/N ratio without this constraint. If the PSFs of 
mirrors were as small as specified in the original design, 
we could also avoid this procedure.) The \ 2 was defined 
as 

(Dt — (Distance)) 21 



E 

T=2,3 



W4 



Distance 




where Wt is the width seen from the intersection point, 
and Dt is the distance between the image center of 
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Fig. 8. — 9 2 distribution for the standard analysis of Crab 
Nebula data. The points with error bars are from the ON-source 
region and the solid histogram is the background region. 



gravity and the intersection for each telescope (Fig. 0), 
(Distance) is the mean distance obtained by Monte- 
Carlo simulations for gamma-rays, and. o Distance 

is its 

standard deviation. The improved 6 2 distribution is 
shown as the solid line in Fig. |SJ The number of events 
with 2 < 0.05 degree 2 increases by a factor of 1.8. Note 
that the angular resolution was estimated to be 0.23 de- 
gree and 0.23 2 is roughly 0.05 degree 2 . 

We first used the conventional ( "square cuts" ) method 
as follows; we accepted events with width <0.2 (for T2) 
and 0.15 (for T3), and length<0.3 (for T2) and 0.25 
(for T3), respectively, where all values are in units of 
degrees. The cut values differ for T2 and T3 due to 
the PSF differences. As the 9 resolution is 0.23 de- 
gree, we took six background points on the 0.5-degree- 
radius circle from the pointing direction. The result- 
ing 9 2 distributions for the ON-source points (the points 
with error bars) and the OFF-source points (the solid 
histogram) are shown in Fig. EI The background normal- 
ization factor is 1/6, equal to the inverse of the number 
of background points in these wobble mode observations. 
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Fig. 10. — Distributions of opening angle and separation angle of 
two images. These were obtained also via the subtraction described 
above for the Crab data. 



Fig. 9. — Hillas moment distributions: upper-left (T2 widths), 
upper-right (T2 lengths), lower- left (T3 widths), and lower-right 
(T3 lengths). The points with error bars were obtained after the 
subtraction of background events. The solid histograms are the 
Monte-Carlo predictions where the total number of entries were 
normalized to those of the observations. 



The number of excess events with 9 2 < 0.05 degree 2 is 
(2207—1990=) 216, where the numbers in parentheses 
are the on-source and estimated b ackground count s, re- 
spectively. This corresponds to a iLi &; Mai ljl983j) sig- 
nificance of 4.4(7. The number of eve nts predicted by 
our M onte-Carlo simulation, using the lAharonian et alJ 
( 2000) flux, an E~ 2 - 5 gamma-ray energy spectrum, and 
minimum and maximum gamma-ray energies of 500 GeV 
and 20 TeV, was 195 events. (For comparison, we also un- 
dertook "mono", i.e., single telescope, analyses using the 
regular CANGAROO-II procedures. Excesses for both 
telescopes were found using both by square cut and Like- 
lihood cut analyses. The statistical significances of these 
excesses were at the 3<r level, confirming the power of the 
stereo technique.) 

Using the events with 9 2 < 0.05 degree 2 , we can derive 
the Hillas moment distributions after background sub- 
tractions. This provides a check of how well the Monte- 
Carlo simulations agree with the real gamma-ray data. 
These are shown in Fig. [5J Our Monte-Carlo simula- 
tions are consistent with the data within statistical un- 
certainty. Note that the width distribution is sensitive 
to the PSF of the mirror system and that that of T3 is 
better than T2. 

After this standard analysis, we investigated two bias- 
free analyses . The first is the Lik elihood method in- 
troduced by lEnomoto et al.1 l)2002|) . In the standard 
"square cuts" , there are four parameters on which cuts 
can be made, and in the absence of a strong gamma-ray 
source or detailed simulations accurately incorporating 
the characteristics of the telescope hardware, some free- 
dom in choosing the exact cuts is a vailable, (see, e.g., 
the discussion in Dazclcy et al. 2001). In the Likelihood 
method, we proceed by making probability density func- 
tions (PDFs) from the histograms of width and length, 
for the two telescopes, and some other parameters such 
as the opening angle and the distance between two im- 
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Fig. 11. — Likelihood Ratio (LR) distributions. The points with 
error bars are from Crab Nebula gamma-rays, the solid histogram 
from the gamma-ray Monte-Carlo, and the dashed histogram from 
all events. The latter two distributions are normalized to the total 
number of the observed excess. 



ages' center of gravity (image separation) (Fig. HOj l. In 
order to treat the energy dependence of these parame- 
ters, we used two-dimensional histograms: the PDFs are 
therefore 2-D functions. For the gamma-ray sample, we 
used the data from Monte-Carlo simulations, and for the 
cosmic ray background sample, we used real observation 
data, since the "contamination" by gamma-rays is much 
less than 1%. Thus for each event, we can derive a prob- 
ability (L) for the event being initiated by a gamma-ray 
or a cosmic ray, where normalization ambiguities still re- 
main. The Likelihood ratio LR is defined as 



LR = 



L(gamma — ray) 



L(background) + L(gamma — ray) 



The distributions for the Crab Nebula data are shown in 
Fig. ^2 The observational data are in reasonable agree- 
ment with the Monte-Carlo predictions. Observationally, 
the optimal cut was found to be LR >0.9. The number 
of excess events is (390—289=) 101 (a Li and Ma sig- 
nificance of 5.2 cr) where the Monte-Carlo expectation is 
121 events. The determination of the cut position is, 
however, still somewhat subjective. 
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Fig. 12. — Statistical significance of the Crab excess versus LR 
cut position (the solid line). The dashed curve is the figure of merit 
calculated as described in the text. The normalization is arbitrary. 
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Fig. 13. — 8 2 distribution for Crab Nebula data obtained using 
the LR > 0.5 cut. 



A simple figure of merit — the number of the ac- 
cepted events in the Monte-Carlo simulation divided by 
the square-root of those in the observation — is plotted 
as a function of LR in Fig. El The Monte-Carlo predic- 
tion is a smoothly increasing function of the cut position 
(the dashed line). On the other hand, the statistical sig- 
nificance of the observed excess is shown by the solid 
line. The two show good agreement over a wide range; 
however, some discrepancy is apparent in the extremely 
tight cut positions. For example, a LR cut at 0.5 gives 
an excess (Fig. Q3|, of (1049-895=) 154 events (4.6tr) 
where the expectation is 182 events. Further fine-tuning 
of the Monte-Carlo code is necessary in the future and 
such tight cuts will for now be avoided. The optimal cut 
position of 0.9 is located on the edge of this tight cut 
region. 



We now investigate an alternative approach for the 
comparison of obser vational and M onte-Carlo data: the 
Fisher Discriminant l)Fisherlll936|) . When we use multi- 
parameters: 

P = (Vector of Hillas moments for two telescopes) 
and assume that a linear combination of 
F = a ■ P, 

provides the best separation between signal and back- 
ground, then the set of linear coefficients (a) should be 
uniquely determined as 



a 



Hsig — HBG 



E m 



E 



BG 



where jj, is a vector of the mean value of P for each 
sample, i.e., = (P), and E is their error matrix, i.e., 
= {PP T ) - (P)(P T ). The values of fl sig , fl BG , E sig , and 
Ebg can be calculated from the Monte-Carlo gamma- 
ray events for signal and observational data for back- 
ground. Our purpose is to separate "sharp images" from 
"smeared ones" . This method is regularly used in high- 
energy physics experiments, s uch as the - B-fac tory's "Su- 
per Fox- Wolfram moment" ijAbe et alJ 12001). in order 
to separate spherical events from jet-like events. Here, 
the assumption strongly depends on which linear com- 
bination is best. We must, therefore, select parameters 
which are similar to each other. Width and length are 
both second order cumulative moments of shower im- 
ages and thus a linear combination is a reasonable as- 
sumption. We used four image parameters: T2-width, 
T3-width, T2-length, and T3-length. Here, the energy 
dependence of width and length for each telescope were 
corrected using the energy estimated from the summa- 
tion of ADC values of hit pixels for that telescope using 
the Monte-Carlo expectations. This correction function 
was a second-order polynomial and we carry out an offset 
correction using it, i.e., 



P = Pr, 



a 



ailog(£2 ADC) - a 2 log(Y J ADCf 



where d^s were determined from the two-dimensional 
plots of raw "Hillas moments" versus log(^2 ADC) us- 
ing the Monte-Carlo gamma-ray simulations and P raw is 
a vector made of the raw values of "Hillas moments", 
respectively. The corrected width and length are dis- 



tributed around zero and the 
cated exactly at zero, 



(Psig) is automatically lo- 



(P sig ) = 

resulting in the Monte-Carlo expectation of F location 
being zero: 



{F„ 



0. 



This removes cut-selection bias: a cut at zero ensures 
50% acceptance (for gamma-rays) automatically, a typi- 
cal value adopted in IACT analyses. From the 8 2 Crab 
Nebula excess we obtain the following distribution of 
Fisher Discriminant, F, as shown in Fig. The be- 
havior of the distributions for gamma-ray signals, back- 
ground events, and Monte Carlo expectations agree rea- 
sonably well and are similar or a little bit better than the 
result obtained for the Likelihood method. This figure 
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Fig. 14. — Distributions of the Fisher Discriminant, F. The 
points with error bars are obtained from the observation, the solid 
histogram from the gamma-ray Monte-Carlo, and the dashed from 
all events representing background behavior. 



Fig. 16.— The final 6 2 plot for the Crab Nebula, using the 
Fisher Discriminant analysis with the cut at zero. 
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Fig. 15. — Statistical significance of the Crab Nebula excess 
versus F cut position (the solid line). The dashed curve is the figure 
of merit calculated as described in the text. The normalization is 
arbitrary. 



should be compared with Fig. 2 in the recently published 
H.E.S.S. pap er for the distribution of their mean reduced 
scaled width ijAharonian et aUl2005b|) . The H.E.S.S. dis- 
crimination is better; however, improved results can be 
expected for CANGAROO-III if the quality of mirrors 
were improved to the levels of H.E.S.S. or MAGIC. 

We consider the figure of merit for this method in 
Fig. 1151 This time, the agreement is very good in 
all regions. Although the cut at F > —0.16 showed 
the best statistical significance, 5.8er, with an excess of 
(1177-974=) 203 events, where the Monte-Carlo expec- 
tation is 162 events. A conservative choice of cut point 
is zero which, as noted above, is exactly the mean ex- 
pected position for a gamma-ray acceptance of 50%. In 



the following we use the Fisher Discriminant analysis as 
the default, with a cut position at zero, thus removing all 
subjective or potentially biased elements from the anal- 
ysis. The excess is (744—604=) 140 events, where the 
expectation is 110 events. The final 9 2 distribution is 
shown in Fig. 1161 

In this analysis, we estimated the gamma-ray energy 
using the average of the T2 and T3 summations of ADC 
values, which probably provides the optimal energy res- 
olution for this analysis. The usage of the intersection 
constraint, due to finite point spread function of the in- 
dividual mirror segments, does not allow a core-distance 
correction of energy. The overall energy resolution was 
considered to be 30% from the simulation study. We 
show the differential flux for the Crab Nebula in Fig. 1171 
The points with error bars are our data where the errors 
are statistical. The systema tic error is considered t o be 
15% at t his stage. HEGRA IjAharonian et alJl2000j) and 
Whipple (Hill as et alJll99q) results are also shown. The 
results agree within the statistical and systematic errors. 
The deviation in the CANGAROO flux at higher energies 
may be due to saturation effects, which are not yet fully 
implemented in our Monte-Carlo code; however, this in 
not an important consideration for Vela observations, 
which are the main focus of this paper. We leave further 
consideration of this for the future, noting that while the 
analysis template defined here for CANGAROO-III is, 
at present, most suitable, it may be able to be further 
optimized with refined Monte-Carlo simulations. 

We show the 2D profile of gamma-ray images in Fig. 1181 
The angular resolution is estimated to be 0.23 degrees 
(ler). Here, only high energy events [E >10.6TeV) are 
plotted in order to improve the S /N ratio and to remove 
uncertainties in the background subtraction. The back- 
ground level at each position was estimated using the 
events with values of F below (Fbg) which are consid- 
ered to be mostly background protons. Here, however, 
we need to note that the subtraction template for the 
background was made using events with low F sample. 
The Fisher Discriminant uses a linear combination of the 
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Fig. 17. — Differential gamma-ray flux from the Crab Neb- 
ula as a function of energ y. The red line is the HEGRA re- 
sult lAh aronia n et alJ I200Q") and the green is the Whipple result 
IHillas et alJH998l) . 



Fig. 19. — The 9 2 distribution for Vela Pulsar observations 
resulting from application of the analysis template described in 
§3.3. 
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Fig. 18. — Excess number plot in the equatorial frame. The 
unit is events per 0.1 X 0.1 degree 2 in the FOV coordinates. The 
cross indicates the center of the Crab Nebula. The resolution is 
considered to be 0.23 degrees in declination in this frame. 



energy-corrected width and length for T2 and T3, respec- 
tively. The subtraction sample for morphological studies 
may tend to be overwhelmed by the larger sized events 
which have poorer angular resolutions. The important 
thing is that the background sample tends to be flatter 
than the gamma-ray events. Therefore, we only show a 
very restricted area in this figure. An improved analysis, 
required for diffuse gamma-ray detection, which will be 
introduced later. 

As a result of the considerations above, the analysis 
template is: 

• the Fisher Discriminant is adopted. 

• the cut position on F is exactly at zero. 



• others parameters, such as elevation cut and 
shower-rate cut, are determined as appropriate for 
the source under investigation (depending on its 
declination and galactic coordinates). 

4. OBSERVATIONS OF THE VELA PULSAR 

The Vela Pulsar was observed between 2004 January 
17 and February 25. In total, the preselected data cor- 
respond to an analyzable period of 1311 min. after the 
elevation and cloud cuts, where the minimum elevation 
angle was set at 60 degrees and the shower rate at 9 Hz. 
The mean elevation angle was 70.9 degrees, correspond- 
ing to an energy threshold of 600 GeV. The observations 
were carried out using the same wobble mode as for the 
Crab Nebula observations. In this period, T2 and T3 
were in operation and we analyzed the stereo data from 
these two telescopes. As the Vela Pulsar is at a declina- 
tion of —45°, the relative orientation of the two telescopes 
does not present any problems. 

We used the optimized analysis procedure described in 
the previous section and so there are no a posteriori tri- 
als to consider in the interpretation of results — it is a 
"blind analysis" . The resulting 9 2 distribution is shown 
in Fig. ["""jl After the background subtraction, we ob- 
tained an excess of (677—722=) — 45±29 events. Monte 
Carlo simulations, with minimum and maximum gamma- 
ray energies of 100 GeV and 20 TeV, respectively, and 
an energy spectrum proportional to E~ 2 ' 5 , predict 394 
events would be detected for a 100% Crab level gamma- 
ray source. The CANGAROO-I flux was 60% of the 
Crab Nebula flux. The F distribution, obtained using 
subtraction, is shown in Fig. 12*01 

There is no excess of gamma-ray-like events around 
the predicted region, offset by 0.13 degrees from the Vela 
Pulsar. Thus observations with significantly improved 
instrumentation and a robust analysis procedure do not 
support the previous claim for TeV gamma-ray emission 
from this region. 

The 2a upper limits are shown in Fig. 1211 together 



10 



Enomoto et al. 



TABLE 1 

2-cr UPPER LIMITS TO THE INTEGRATED FLUXES. 



F(> E) at F(> E) at 

Energy, E offset position pulsar position 
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(cm 2 s 1 
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Fig. 20. — F distributions. The points with error bars are 
obtained from the observations, the solid curve from the gamma- 
ray Monte-Carlo simulations, and the dashed histogram from all 
observed events to represent the cosmic ray background behavior. 
The Monte-Carlo histogram is normalized to a 100% Crab Nebula 
flux level. 
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Fig. 21. — 2a upper limits f or the gamma-ray flux. C-I repre- 
sents the CANGAROO-I result lYoshikoshi et aljll997t) and C-III 
indicates the CANGAROOTII result reported he re. Also shown 
are Durham IChadwick et alJ 120001). BIGRAT ID azelevl I [19991) 
and H .E.S.S. preliminary results fKhelifi et aU20 05) iKhcl ifi et alJ 
20051). 



with the results of other observations. The upper limits 
in the figure are a factor of 5 below the CANGAROO-1 
result. This analysis has used the point offset by 0.13 
degrees to the south-east of the Vela Pulsar, i.e., offset 
is (ARA,Adec)=(0.14°,-0.1°) from the pulsar position. 
This position which was the maximum of the excess de- 
tected with the CANGAROOT telescope: an analysis 
at the position of the Vela Pulsar position yields similar 
upper limits. These are summarized in Tabled 

5. VELA PULSAR WIND NEBULA 



Recently, extended TeV gamma-ray emission coinci- 
dent with the Vela Pulsar Wind Nebula (PWN) wa s 
reported by the H.E.S.S. group (|Khelifi et a.U l2005h|L 
Their preliminary report claimed that the center of the 
emission is (RA, dec) = (8 h 35 m , -45°36') and that the 
flux within a 0.6 degree radius from this position is 50% 
of the Crab Nebula at 1 TeV. They also noted that no 
gamma-rays were detected from the Vela Pulsar, and 
placed a tight upper limit on the pulsar flux at 0.26 TeV 
ijKnelifi et a,l]l2005lhlL 

Thus far, we have focused on a point source analy- 
sis based on the wobble mode observations. The peak 
PWN source position coincides with the one of two wob- 
ble pointing direction in the coordinate of the field of 
view of the camera. We, therefore, can not carry out 
background estimation using the usual wobble method. 
In this section we undertake an optimized analysis for 
an extended source. Another difficulty is that we don't 
have sufficient statistics for OFF source regions, and so 
the background subtractions should be carried out us- 
ing the ON source data runs. Gamma-ray-like events 
can be extracted by fitting position-by-position F dis- 
tributions under the assumption that gamma-rays obey 
the Monte-Carlo predictions, the proton background fol- 
lows the average F distribution of all directions, and the 
total distribution is a linear combination of those two. 
The separation between those two distributions is likely 
to be worse at lower energies due to the smaller image 
sizes. These distributions are plotted in Figs for var- 
ious energy ranges. From the upper to lower panels, the 
central energies of the energy bins are 540, 780, 1200, 
2700, 3300, and 8800 GeV, respectively. The blank his- 
tograms are gamma-rays and the hatched, protons. As 
shown in figure, the separation begins at 780 GeV and 
becomes significant at 1200 GeV. Therefore, we first an- 
alyzed events with energy greater than 1200 GeV. 

Then we checked the directional dependence of the 
F distributions. The field of view was segmented into 
0.2x0.2 degree 2 regions and each F distribution was 
compared with the average. The reduced x 2 distribu- 
tion is shown by the histogram in Fig. 1231 The curve 
shown is the predicted distribution for the 23 degrees 
of freedom. There is good agreement, i.e., there is no 
significant directional dependence. 

This method was then checked with Crab Nebula data. 
The F distributions at various 9 2 slices were taken at en- 
ergies greater than 5.7 TeV. The background F distribu- 
tion was obtained in the higher 8 2 region, 0.1-0.3 degree 2 . 
The only fitting parameter is the percentage of gamma- 
ray-like events relative to the total events. The result is 
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Fig. 22. — 7* 1 distributions in various energy regions. The blank 
histograms are for Monte Carlo simulations of gamma-rays and the 
hatched histograms for background regions of real data. From the 
top panel to the bottom, the central gamma-ray energies are 540, 
780, 1200, 2700, 3300, and 8800 GeV, respectively. The histogram 
entries are normalized to one event in the shown regions. 



Fig. 24. — The 9 2 plot obtained by the fitting method de- 
scribed in the text for Crab Nebula data. The vertical scale is the 
percentage of gamma-ray— like events relative to all events. 
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Fig. 23. — Reduced \ distribution for various arrival di- 
rections (the histogram). The pixel size for each direction is 
0.2 X 0.2 degree 2 . The curve is the expected distribution for the 
number of degrees of freedom, 23. 
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Fig. 25. — The wide range 6 plot for the Vela PWN region, 
where 9 was calculated from (RA, dec) = (8 h 35 m , -45° 36'), i.e., 
the peak of the emission detected by H.E.S.S. The gamma-ray- 
like events were extracted by the fitting procedure described in the 
text. 



shown in Fig. [21 The statistical significance of the peak 
is 4.0cr, while the ordinary wobble analysis gave a 3.6c 
excess. 

Having demonstrated the validity of this method, we 
carried out an analysis of the Vela PWN region. The 
H.E.S.S. group detected a gamma-ray excess extended 
over a 0.6 degree radius from the center of the emission 
[(RA, dec) = (8 h 35 m , -45°36')]. In our case, the angular 
resolution (0.23 degree) is significantly worse than that 
of H.E.S.S. We, therefore, chose the background region 
to be more than 0.8 degree from the center. The result 
of fitting is shown in Fig. [23 An excess was observed at 



2 < 0.6 degree 2 around the center of Vela PWN region. 
The excess radius is marginally consistent with H.E.S.S. 
considering our angular resolution. The total number of 
gamma-ray-like events is 561±114. 

The differential flux was obtained and is shown in 
Fig. [213 Although the statistics are poor, the spectrum 
looks hard, consistent with the preliminary H.E.S.S. re- 
sults, and the fluxes are in general agreement. This ex- 
cess is bel ow the 5 a ley el generally required for a firm 
detection ijWeekes |l2005j) and is spatially extended over 
a significant portion of our field of view, and thus a high 
resolution (0.1 degree) morphology map is not justified. 
We are, however, able to offer supporting evidence of the 
H.E.S.S. result and, in this light, are planning to observe 
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Fig. 26. — Gamma-ray flux in the Vela PWN region. 



stereoscopic imaging Cherenkov telescopes. At that time, 
two telescopes (T2 and T3) were in operation and events 
coincident to the two telescopes were analyzed. Cali- 
bration was performed using muon rings and the per- 
formance of the telescopes confirmed with observations 
of the Crab Nebula. The estimated energy threshold 
for this analysis was 600 GeV. The use of the Fisher 
Discriminant has been introduced and a template for 
CANGAROO-III analysis presented. No significant ex- 
cess of events was found from the Vela Pulsar direc- 
tion or from the peak of the emission detected with the 
CANGAROO-I telescope. The upper limits obtained are 
a factor of five less than the CANGAROO-I fluxes (as- 
suming an _E~ 2 - 5 spectrum). In addition, we have con- 
firmed, at the 4cr level the gamma-ray emission recently 
reported from the Vela Pulsar Wind Nebula. The TeV 
emission from the PWN peaks ~0.5 degree south of the 
pulsar position and has a ^0.6 degree extension, consis- 
tent with the H.E.S.S. report. A detailed morphological 
study of this source will require more observations. 



the PWN region in more detail next year. 

6. CONCLUSION 

We have observed the Vela Pulsar region from 2004 
January 17 to February 25 with the CANGAROO-III 
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APPENDIX 
CRAB OBSERVATION 

Observations of the Crab Nebula were carried out over the period 2003 December 18 to 28. The observations were 
made using the so-called "wobble mode" in which the pointing position of eac h telescope was al ternated between 
±0.5 degree in declination from the center of the Crab Nebula every 20 minutes ijDaum et al.lll997j) . Two telescopes, 
T2 and T3, were operational at this time. T3 is located 100 m to the south-south-east of T2, which is not ideal for 
observations of northern sources; however sufficient data was recorded for a useful analysis to be made. The trigger 
rates of the individual telescopes were at most 80 Hz and this was reduced to 8 Hz by requiring the above coincidence. 
The total observation time was 1122 min. Next we required both telescopes to have clusters of PMTs with five adjacent 
hits above a 5p.e. hit threshold, which reduced the event rate to 6 Hz. Looking at the time dependence of this event 
rate, we can remove data taken in cloudy con ditions. This procedure is exactly the same as the "cloud cut" used in 
CANGAROO-II analysis ijEnomoto et aljf2002|) . Only data taken at elevation angles greater than 30 degrees were 
used, resulting in a total of 890 minutes data, with a mean elevation angle of 35 degrees. 
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